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1. Introduction 

In 1918 Emmy Noether wrote the seminal paper "Invariante Variationsprobleme" [I], 
where she showed that for differential systems derived from a variational principle, 
conservation laws could be obtained from Lie group actions which left the functional 
invariant. 

Recently in [2] , it was proved that for Lagrangians that are invariant under some 
Lie symmetry group, Noether's conservation laws can be written in terms of a moving 
frame and vectors of invariants. Furthermore, in .2] the authors showed that for one- 
dimensional Lagrangians that are invariant under a semisimple Lie symmetry group, 
the new format for the conservation laws could reduce considerably the calculations 
needed to solve for the extremal curves. In particular, a classification is given for 
variational problems which are invariant under the three inequivalent SL(2, C) actions 
on the plane (classified by Lie [3]). In this paper we use the new structure of Noether's 
conservation laws presented in |2] to simplify variational problems that are invariant 
under SE{2) and SE(3); these groups are not semisimple. 

In Section [2j we will briefly give an overview on moving frames, on differential 
invariants of a group action and on invariant calculus of variations. Throughout 
Section [2] we will use the group action of SE(2) on the plane as our pedagogical 
example. 

In Section[3J we show in some detail how to compute the new version of Noether's 
conservation laws presented in Theorem 12.91 for one-dimensional variational problems 



Moving frames and conservation laws for Euclidean invariant Lagrangians 



2 



that are invariant under SE(2), and then demonstrate how their invariantized 
Euler-Lagrange equations and Noether's conservation laws can be used to solve the 
integration problem. 

Finally in Section^ we present the simplified solution to the physically important 
one-dimensional variational problems that are left unchanged under the SE(3) group 
action. 

2. Structure of Noether's Conservation Laws 

In this section, we will give a brief overview of concepts regarding moving frames, 
differential invariants of a group action and the invariant calculus of variations needed 
to understand the statements of our result. For more information on these subjects, 
see Fels and Olver [US], Mansfield [6] and Kogan and Olver [7]. We will use the SE(2) 
action on the plane as our pedagogical example. 

2.1. Moving frames and differential invariants of a group action 

Here we are using moving frames as reformulated by Fels and Olver [H[S], adapted to 
the context of differential algebra. 

Let X be the space of independent variables with coordinates x = {xx, x p ) and 
U the space of dependent variables with coordinates u = (u 1 , ...,u q ). We will use a 
multiindex notation to represent the derivatives of u a , e.g. 

g\K\ u a 

II 

K dx^dx? ■■■dxp p ' 

where the tuple K = [kx, k p ), represents a multiindex of differentiation of order 

\K\ = kx + k 2 H + kp. Hence, let M = J n (X x U) be the n-th jet bundle with 

coordinates 

z = (xx, ...,x p ,u x , ...,u q ,u\, ...). 
On this space, the operator d/dxi extends to the total differentiation operator 

a = l K A 

In this paper we are interested in using Noether's conservation laws to find the 
solutions that extremize variational problems which are invariant under SE(2) and 
SE(3). Thus, consider a Lagrangian L(z)dx that is invariant under some symmetry 
group. Let a group G act on the space M as follows 

G x M -> M 

z i— > Iz = g ■ z, 

which satisfies either g ■ [h ■ z) = (gh) ■ z, called a left action, or g ■ (h • z) = (hg) ■ z, 
called a right action. We say a Lagrangian L(z)dx is invariant under some group 
action if 

L(z)dx = L(z)dx 

for all g G G. 

Consider a Lie group G acting smootlhy on M such that the action is free and 
regular. Then for every z £ M there exists a neighbourhood U of z, as illustrated in 
Figure Q] such that 
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- the group orbits have the dimension of the group G and folliate U ; 

- there is a surface K, C U which intersects the group orbits transversally at a single 
point. This surface is called the cross section; 

- if O(z) represents the group orbit through z, then the group element g <E G taking 
z £ Li to k is unique. 

U 

K 

S\ ^ O(z) 

z 




Figure 1. A local foliation with a transverse cross section 

Under these conditions, we can define a right moving frame as the map p : U — > G 
which sends z£Wto the unique element p(z) 6 G such that 

p(z)-z = k, {k} = o(z)nic. 

The element g € G in Figure [T] corresponds to p(z). 

To obtain the right moving frame, which sends z to k, we must first define the 
cross section K. as the locus of the set of equations tpj(z) = 0, j = 1, r, where r is 
the dimension of G. Normally, the cross section is chosen so as to ease the calculations. 
Then solving the set of equations 

= i/jj(g- z) = 0, j = l,...,r, (1) 

known as the normalization equations, for the r group parameters describing G yields 
the right moving frame in parametric form. Hence, the frame obtained satisfies 

i>j{p{z) ■ z) = 0, j = l,...,r. 

By the implicit function theorem, a unique solution of ([T]) provides an equivariant 
map, i.e. for a left action 

p(z) = p{z)g~ 1 

and for a right action 

p(z) = ,9~V(z). 

Example 2.1. Consider the SE(2) group acting on curves in the (x,y(x)) -plane as 
follows, 

( x \ / x \ / cos 9 sin 9 \ / a; — a 
V 2/ / ^ V V ) ~ \ ~ s ' m0 cos6 J \ y-b 

where 9, a and b are constants that parametrize the group action. Here we are using 

the inverse action because it simplifies the calculations. 
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There is an induced action on the derivatives yx, where K is the index of 
differentiation with respect to x, called the prolonged action. The induced action 
on y x is defined to be 

„ dy dy / dx 

y x = g-y x = j~ = 

by the chain rule, so the action of flp on y x is 
- sin 9 + y x cos 9 



!J.r 



cos 9 + y x sin 9 



Similarly, 



_ d 2 y _ 1 d / dy\ _ y xx 
Vxx di.r;- ~ dx/dxdx \dxj~ (cos9 + y x sm9) 3 - { ' 

If we consider M to be the space with coordinates (x,y,y x ,y xx , ...), then the action is 
locally free near the identity of SE(2). Thus, taking the normalization equations to be 
x — 0, y — and y x = 0, we obtain 

a — x, b = y, and 9 = arctanj/ x (4) 

as the frame in parametric form. 

Remark 2.2. In this paper we will consider all independent variables to be invariant. 
If these are not invariant, then we can reparametrize and set the original independent 
variables as depending on the new invariant parameters. 

Theorem 2.3. Let p(z) be a right moving frame. Then the quantity I(z) — p{z) ■ z 
is an invariant of the group action (see J^j). 

Consider z — (z%, z m ) G M and let the normalization equations Zi — Ci for 
i = 1, r, where r is the dimension of the group G, then 

p(z) ■ z = (ci, ...,c r ,I(z r+1 ), ...,I(z m )), 

where 

I(zi) = g ■ z\ g=p{z) , l = r + l,...,m. 

Example 2.1 (cont.) Evaluating yfff, given by (0) at the frame ^ yields 

1 Vxx 

J/xx|(a=a;,6=a,e=arctanj/ x ) — 7^ ~ y2")3/2 ' 

the Euclidean curvature. So evaluating y~K at the frame yields a differential 
invariant. 

Definition 2.4. For any prolonged action in the jet space J n (Xx U), the invariantized 
jet coordinates are denoted as 

Ji = I(xi) = x,\ g=p(z) , 1% = I{u a K ) = u%\ g=p(z) . (5) 

These are also known as the normalized differential invariants. 

Example 2.5. Consider the group action of SE{2) as in Example \2.1\ Since x is 
not invariant we reparametrize (x,y(x)) as (x(s),y(s)), where s is invariant and let 
g G SE(2) act on (x(s),y(s)) as in Example \2.1\ Solving the normalization equations 
x = 0, y = and y s — 0, we obtain the frame 

a — x, b = y, 9 = arctan [ — ) . (6) 
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We have then 

9- z \g= P {z) = (s,x,y,x s ,y s ,yZ)\ g=p ^ 
= {I{ S )J*J\I?JIJ1 1 ) 

= L o, o, yffiTyl, o, x ' y ''7"' x '' ) ■ (7) 

V V X s +Vs J 

The second, third and fifth components of correspond to the normalization 
equations x — 0, y — and y s — respectively. The fourth and sixth components, If 
and I\ x respectively, are the lowest order differential invariants and all higher order 
invariants can be obtained in terms of them and their derivatives. 

Theorem 2.6. (Replacement Theorem J5j/j If f(z) is an invariant, then 

f(z) = f(I(z)). 

This theorem allows one to find the ij£ in terms of historically well-known 
invariants without having to solve for the frame. 

Example 12.51 (cont.) We know that SE(2) preserves |x s |, thus applying the 
Replacement Theorem we obtain 



which yields that \J x 2 s + y 2 = If up to a sign. Next we know that the Euclidean 
curvature n is also invariant under SE{2), and we obtain 

_ x s y ss y s x ss ^ ^1^11 iflii ^11 

K - (^ + y 2 ) 3/2 ~ ((W+IWW ~ WF 

which gives us I\ x in terms of k and |x s |. 

Since we are considering all independent variables to be invariant, all total 
differential operators will also be invariant. Hence, the invariantized differential 
operators 

We know that 
d 

——n a — 11° 

dx n K -u Ki , 

although the same cannot be said about its invariantized version and in general 
T>ilx ^ I^; indeed we have that 

^ilK = lKi+Mh, (8) 

where is known as the correction term. We will not go into its calculation since it 
would take us to far afield, but for more information on correction terms see §4.5 [BJ. In 
any case, software exists to calculate these correction terms [8]. Equation [8] shows that 
the processes of invariantization and differentiation do not commute. Considering two 
generating differential invariants If and /£ and letting JK = LM so that If K = I" M , 
then this implies that 

V K I a j - Mj K = V M It - M£ M . (9) 
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These equations are called syzygies or differential identities. These will play a crucial 
role in the obtention of the invariantized Euler-Lagrange equations and Noether's 
conservation laws. 

Example 12.51 (cont.) // we set x = x(s,t) and y — y(s, t) and take the normalization 
equations as before, we obtain 



x s x t + y s yt ~, 
xt\ g =p(z) = h = — 75f=f=r> yt\g=p(z) 



_ T y _ x sVt ~ VsXt 



Furthermore, since both s and t are invariant, T> s and T>t commute. From Figure 
we can see that there are two ways in which we can obtain If 2 



1 T1 



JX TX TX m 

If Ifi • 

Figure 2. Paths to If 2 



and since both ways must be equal, we get a syzygy between 1% and rj 
is 



V t 7] = VJl - KTjIl 

Similarly, we have a syzygy between 1\ and F\ x and the syzygy is 

Hi 



2nr[D s Il + n s rjl%. 



The syzygy 
(10) 

(11) 



2.2. Invariant calculus of variations 

Kogan and Olver in [7] studied invariant calculus of variations from a geometric point 
of view, we instead do it from a differential algebra point of view. 

Assume Lagrangians to be smooth functions of x, u and finitely many derivatives 
of the u a and denote these as -Sf[u] = J L[u]dx. Furthermore, suppose these are 
invariant under some group action and let the kj, j = 1, ...,7V, denote the generating 
differential invariants of that group action. Also, assume that the action leaves the 
independent variables invariant so that the Lagrangians can be rewritten as J L[/c]dx. 

To obtain the invariantized Euler-Lagrange equations we proceed in a similar way 
as for finding the Euler-Lagrange equations in the original variables (x, u). 

Recall that if x 4 (x, u) extremizes the functional J£f[u], then for a small 
perturbation of u, 

0=4^- Jz^u + ev] 



ds 



Q 

E 

a=l 



E a (L)v a +Y—(^v a 
v ' 4- DxA duf 



dx 
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after differentiation under the integral sign and integration by parts, where 

D \K\ Q 



K 



Dx 



fei 



■ Dx, 



du 



K 



is the Euler operator with respect to the dependent variables u a . The boundary terms 
correspond to Noether's conservation laws and the variation v to the infinitesimals. 

To obtain the invariantized Euler-Lagrange equations, we first introduce a dummy 
invariant independent variable t and set the u a — u a (x,t). The introduction of this 



new independent variable results in q new invariants 1° = g 
syzygies V t K — V.I(u t ) that is 



a 



t \g=p(z) 















1- 


[") 

Uv 



and a set of 



(12) 



where % is a N x q matrix of operators depending only on the T>i, for i = 1, ...,p, the 
Kj, for j = 1, N, and their invariant derivatives. Since all independent variables are 
invariant, we have that all differential operators commute, specifically [T>i,T> t ] = 0, 
for all i — 1, ...,p. 

Remark 2.7. Up to this moment we have represented the independent variables as 
x and the dependent variables as u. Since the examples in this paper only involve 
two independent variables, s and t, we will represent the coordinates of the space 
of dependent variables as x, where the dependent variables will be the usual space 
coordinates. 

For simplicity, consider a one-dimensional Lagrangian L{x, y, y x , y xx , ...)dx with 
a finite number of arguments which is invariant under the SE(2) group action ([2]). 
Such variational problems can be rewritten in terms of the generating invariants of its 
group action, in this case the Euclidean curvature, k, and its derivatives with respect 
to s, the Euclidean arc length. We reparametrize (x,y(x)) as (x(s),y(s)), and to fix 
parametrization as arc length, we introduce r\ = \J x^ + \ 
we consider the invariantized variational problem 



1 as a constraint. Thus, 



[L(k, k s , k ss , ...) - A(s)(?7 - l)]ds, 



(13) 



where A(s) is a Lagrange multiplier. This constraint does not reduce the solution set 
and it will simplify the calculations. Symbolically, we know that 



d 
d7 



JSf [x + ev] 



e=0 



D 
Dt 



[x]. 



Hence, after differentiating (|13[) under the integral sign and integrating by parts we 
obtain 



T> t / [L(k, Kg, k ss , ...) - X(s)(n - l)]ds 



— V t K + — V a V t K + - — V 2 s V tK 

OK OK, OK„ 



A(s)D t ry 



ds 
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m— 1 n— 



/w, 



ds 



— ^ ^ m=l n=0 ^ ^ 



ds, 



where 



ds m ' 



Next we substitute the circled T> t r\ and 2?tre by their respective syzygies 



V s -K 
K e Vl + K 2 



I V 



Ho 



I(*t), (14) 



where we have already set r\ — 1; it makes no difference to the final result to do this 
at this point. Hence, 

(E K (L)%-A( S )Wx)/(x t )+P s (^ ^(-l) n © s " [j— )v™-^V tK )]ds 

m=l n=0 ^ Km ' 

J [ (H* 2 (E K (L)) - Ht(X(s))) /(x t ) + V s ( - A(s)/f + E K (L)P S /| - V s E«(L)q 



V ™~l-n VtK 



ds, 



(15) 



m— 1 n— 

after a second set of integration by parts, and where HI and H 2 are the respective 
adjoint operators of Hi and H 2 - The vector I(x t ) corresponds to the variation v and 
thus the coefficients of I 2 and I 2 represent respectively the Euler-Lagrange equations 

E X {L) = Ks E K (L) + \ s , (16) 
E y (L) = V 2 S E K (L) + k 2 E k (L) + X(s)k. (17) 

We will use E X (L) = to eliminate A(s) instead of E V {L) = 0, as it contains 
derivatives of re of lower order. Since L does not depend on s explicitly, by the result 
in page 220 of [5], re s E K (L) is a total derivative, specifically 



m—1 j— 



then we obtain 



m— 1 

A(s) = -L+£ £(-l) J '^ 

m= 1 j— 



(18) 



where the constant of integration has been absorbed into A(s) (see Remark 7.1.9. 
of [6]). Hence, we are left with one invariantizcd Euler-Lagrange equation in one 
unknown 



dL 



E"(L) = V 2 S E*(L) + k 2 E k (L) -* £-££ (-1)™^' 

\ m=l j=0 

This particular equation agrees with the one appearing in Kogan and Olver [7]. 



(19) 
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2.3. New version of Noether's conservation laws 

If one calculates the conservation laws for one-dimensional variational problems that 
are invariant under SE{2) from Noether's First Theorem (for the formulae of these, 
see Theorem 4.29 of [9]) and then rewrite these in terms of the Euclidean curvature k 
and its derivatives with respect to the Euclidean arc length s, one obtains 

-A(s) - nE K {L) 

-V S E K (L) 1 = c. (20) 
g= P (z)-^ \ E K {L) 




— sin ft 
cos 9 
a cos 9 + b sin 9 



R{g) v(I) 

where p(z)~ 1 is the right moving frame A(s) is the Lagrange multiplier obtained 
in (|T5)) and c the constant vector. 

We note that R(g) is the Adjoint representation of SE(2) on its Lie algebra se(2). 
To see how calculations in Section 12.21 can yield the result in (f2TJ)) , we first show how 
the Adjoint representation is calculated in the context we will need. 

Consider the SE{2) group action 



x \ _ I cost/ — sine/ \ I x \ ( a 
y ) = { sin0 cos9 ){y ) + {b 

with generating infinitesimal vector fields 

d x , d y , -yd x + xd y . 

Let g e SE{2) act on 

v = ad x + (3d y + j(-yd x + xd y ) 

as in (|2ip. where a, (3 and 7 are constants. Hence, 

g ■ v = adx + Pdy + j(~yds + xdy) 



(21) 



= ( a /3 7 ) 



cos v — sin t) 

sin 9 cos 9 

a sin 9 — b cos 9 a cos 9 + b sin 9 




Thus, R(g) is the Adjoint representation of SE(2), denoted as Ad(g). 

Remark 2.8. In Example \2.1\ we used the right action of SE(2) on the plane to 
calculate the right moving frame, as it simplified its calculation. However, to compute 
Ad(p)^ 1 we considered the left action of SE(2) on the plane, avoiding in this way the 
need to calculate the inverse of Ad(p). 

Next recall the boundary terms (| 15[) obtained in the calculation of the 
invariantized Euler-Lagrange equations, 

m — 1 



A( S )/| + E*(L)v s iy - v s E*(L)iy + E E ( jr- ) = c 

1 n V OK r 



where c is a constant. Substituting 2? s /f and 1 n T> t K for all m in the above 
expression by the differential formulae 

T) TV _ TV _ K pc 
T>tK = I\i2 — 2k/ 12 , 

DsDt^ — -^1112 ^12 ^^-^112 ^^s^12i 
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where these were obtained from |8]). and rewriting it as 



/ -A(«) - «E"(L) \ 



/pi rz 
V J 2 J l 



2k 



dL 

s dn ss 



" 3 ^ + 



/ -V S E K (L) 
E«(L)-2 K 2 ^ 



"(«2 







dL 



V 



fc (22) 



/ 



c« 



yields the boundary terms in a form that is linear in the I% K . 

We now let i be a group parameter. If the parameters are (a\,...,a r 
then from Theorem 3 of [2], it is shown that the vectors ( 1% /f 2 



and 



t 



3 ' 



can be written as the product of row j of Ad(p) 1 and the matrix of invariantized 
infinitesimals 

n-(i) = (cj(D) , g = g 

where a represents a dependent variable, a,j a group parameter, and e the identity 
element. The vector of invariants in (|20p equals the sum of the products of the matrices 
of invariantized infinitesimals fl a (I) with the vectors C a , 

I -X(s)-kE k (L) ' 
O x (I)C a: + n v (I)C v = -V S E K {L) 
V E"(L) 

where 



Q x (7) 



a / 1 
6|00 

-K 



a /'0 
fF(7) = 6 I 1 



Vsss 






1 







Noether's conservation laws for one-dimensional Lagrangians invariant under 
SE(2) can be written as 

\ / -A(s) - kE k (L) 



x s 

Vs 

xy s - yx s 



-Vs 



xx s + yy s 




-V S E K (L) 
E»(L) 



(23) 



where 



yl = 1 was used to simplify the conservation laws and 



A(*) = -L- 



m=l j=0 



dL 



^m—j • 



The following theorem generalizes what we have just seen for one-dimensional 
variational problems that are invariant under SE(2); it states that the conservation 
laws from Noether's First Theorem can be written as the divergence of the product of 
a moving frame with vectors of invariants. 

Theorem 2.9. Let J L(kx> «2> ...)dx be invariant under G x M — > M, where M = 
J n (X x U), with generating invariants Kj, for j = 1,...,N, and let 
for i = I,..., p. Let (ai,...,a r ) be coordinates of G near the identity e, and Vi, for 
i = 1, ...,r, the associated infinitesimal vector fields. Furthermore, let Ad(g) be the 
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Adjoint representation of G with respect to these vector fields. For each dependent 
variable, define the matrix of infinitesimals to be 

where 

A _ dzj 
3 <9a, 

9- 

are the infinitesimals of the prolonged group action. Let fl a (I), for a = 1, q, be the 
invariantized version of the above matrices. 

Introduce a dummy invariant variable t to effect the variation and then integration 
by parts yields 

d 



J L( Kl ,K 2 ,...)dx= J [£E a (L)J? + Div(P) 



dx, 



dt 

where this defines the p-tuple P , whose components are of the form 

Pi = ^2 ItKC^Ki ' = 1; -'Pi 

U,.J 

and the vectors Cf = {Cf K ). Let p(z)^ 1 be a right frame with canonical invariants 
I" K = I(u" K ), where K is the index of differentiation with respect to the independent 
variables Xi, for i — 1, ...,p. Then the r conservation laws obtained via Noether's First 
Theorem can be written in the form 

J2^-Ad(p(z))- 1 v i (I)=0, (24) 



where 



v l (i) = yn a (i)c?. 



The proof can be found in [2]. 
Remark 2.10. One can notice that none of the summands of 

vri—l n— ^ 

in the boundary terms of H15\) show up in the conservation laws i23\) . in other words 
all boundary terms coming from the first set of integration by parts have disappeared. 
This is no coincidence, it is due to the conflation of t with each group parameter; the 
proof of this can be found in )1 0\j . 



In this paper we are interested in showing how the structure of Noether's 
conservation laws can be used to solve the integration problem for variational problems 
that are invariant under SE{2) and SE(2>). In the present section we have computed 
Noether's conservation laws for one-dimensional Lagrangians that are invariant under 
SE(2). In the next section, we will see how these can be used to solve the extremization 
problems. 
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3. Invariant Lagrangians under SE(2) 

Recall that the invariantized Euler-Lagrange equation for a one-dimensional 
Lagrangian invariant under SE(2) is 

E"(L) = V^{L) + k 2 E«(L) -Ji-EE H)'^ (|r) I ' ( 25 ) 

\ m—1 j—0 \ m/ j 

and its associated conservation laws are 

x s -y s \ / -\{s)-kE k {L) 

Vs x s -£> s E re (L) 

xy s -yx s xx s +yy s 1 / V E K (L) 



where 

m — 1 



a( s) = -l+eE(- 1 )^ / /?I 



m=l j=0 

Multiplying both sides of Ad(p(z))~ 1 v(I) = c by Ad(p(zj) yields the following 
system of equations 

- A(s) - kE k (L) = x sCl + y s c 2 , (26) 

-V S E K {L) =x s C2-y sCl , (27) 

E K (L)=y Cl -xc 2 + c 3 . (28) 

Also, we obtain a first integral of the Euler-Lagrange equation as follows. Define 

/ 1 
B= 1 

V o o o 

which satisfies the following equality 

B> = Ad{p)- T EAd{ P y 1 . 
The first integral of the Euler-Lagrange equation is then 

v T (I)Bv(I) = c T Bc 

i.e. 

(A(s) + kE k {L)) 2 + {V S E K (L)) 2 = cj + c\. (29) 

Once k is known, one can see that the integration problem has a straightforward 
solution. Hence, integrating both sides of Equation (121)1) with respect to s yields 



XC\ 



yc 2 = j {-\(s)~KE K (L)}ds, (30) 

and thus solving the two linear equations (|28l) and (I3D1) with respect to x and y one 
obtains 

x ( s ) = ^T^2 ( c i / h A ( s ) - kE k (L)] As - c 2 E K (L) + c 2 c 3 ) , (31) 
y(s) = -^-^ (c 2 f [-X(s) - kE k (L)] ds + Cl E K (L) + . (32) 

C l + C 2 V J C l / 
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One can easily see that Equation (|27[) is immediately satisfied as it is the derivative 
with respect to s of P5]), 

If we consider the Lagrangian with L = 1 and plug it in Equations (j2"ij)) . ipS]). 
(|3ip and (|3"!?1) . then the solution which minimizes the arc length is the equation of a 
line, as one would expect. Another famous Lagrangian to consider is 

J n 2 ds. 

Using the equations above, we obtain, as Euler himself did [TT], that the curvature of 
the minimizing curve satisfies 
1 3 _ 

or the first integral of the Euler-Lagrange equation 
Ak 2 s + n i = c? + ef, 

which is solved by an elliptic function. Solutions are known as Euler 's elastica. For a 
good historical report see [12] . 



4. Invariant Lagrangians under SE(3) 

In the previous section we showed that for an invariant Lagrangian under SE{2), the 
invariantized Euler-Lagrange equation and its associated conservation laws could be 
used to solve the extremising problem. In this section we will proceed analogously and 
present the solution to one-dimensional variational problems that are invariant under 
the following SE{2>) group action on the (x(s), y(s), z(s))-space parametrized by the 
Euclidean arc length s 

x^ x = FT^x-a), (33) 

where R _1 represents the rotation in the three-dimensional space 

cos /3 cos 7 cos f3 sin 7 sin f3 

- sin a sin /? cos 7 — cos a sin 7 — sin a sin f3 sin 7 + cos a cos 7 sin a cos /3 

- cos a sin (3 cos 7 + sin a sin 7 — cos a sin /3 sin 7 — sin a cos 7 cos a cos ft 

and a = ( a b c ) T the translation vector with a, /?, 7, a, b and c as the constants 
that parametrize the group action. 

To solve SE(3) invariant variational problems, we need to find the element 
g E SE(3) that sends the tangent to the curve to the x-axis, the normal to the 
curve to the y-axis and the point (x 7 y, z) to the origin, in other words which sends 
z = (x,y, z,y Sl z s , z ss ) to the cross section (0,0,0,0,0,0). For that we solve the 
normalization equations x — 0, y = 0, z = 0, y s ^ 0, z s = and z^ s — 0, and 
thus obtain the right moving frame 



a — x, b = y, c — z, a = tan 



1 / y,(v s z„-z,y S s)-x,( 



z^x * * — x s z s 



^Jx 2 3 +y^+zl(x s y ss -y 3 

). 7 = tan-i(^). 



(34) 



Consider a one-dimensional variational problem Jz? [x] that is invariant under the 
group action (1331) . To obtain the invariantized Euler-Lagrange equations, we first 
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rewrite jSf [x] in terms of the generating invariants of the group action, which are the 
Euclidean curvature 

II x s x x ss || 



and torsion 

Xsss ' (Xs X X ss ) 

1 1 x s x x ss 1 1 2 

and their derivatives with respect to s. 

Since s represents the Euclidean arc length, the constraint n = yjx 2 + y 2 + z 2 = 1 
must be introduced into the variational problem in order to fix parametrization. 
Hence, the resulting invariantized functional is 

[L(k, t, k s , t s , k ss ,t ss , ...) - A(s)(r? - 1)] ds, (35) 

where A(s) is a Lagrange multiplier. As for SE(2), this will not reduce the solution 
set and will simplify the computation of the conservation laws. 

Next we introduce a dummy invariant independent variable t and set x = x(s, t) 
to effect variation. The introduction of a new independent variable results in three 
new invariants for a = a;, y, z, and a set of syzygies 

In 

H\ k | . C-!(,i 
where the matrix of operators % is 

/ V s — k \ 

Ks K 2 -T 2 + V 2 S -T S - 2tV s 

\ts+2tV s T>s(^) + (^-^)v s + ^V 2 Vs(-^) + (k-^)Vs-^V 2 s + ±VI)' 

where we have already set 77 = 1. 

As in Section 12.21 we differentiate (1551 with respect to t and then integrate by 
parts twice to obtain the invariantized Euler-Lagrange equations 

E X (L) = k s E K (L) + t s E t (L) - V s (2tE t (L)) + A s , 

E*(L) = 2? 2 E K (L) + —VlU{V) +(-- ^) V S E T (L) + (k 2 - r 2 )E"(L) 

+ 2thE t (L) + A(s)k, 
E*(L) = - W(L) + ^P 2 E^(L) + + T - - 2 -4 - k) V s F{L) 




n s E T {L) + 2rV s E K (L) + r s E K (L), 
and the boundary terms 

(2rE T (L) - A) 7f + (—Z T {L) - —VsE^L) - V a E K {L)) J* + f E K (L) + — E T (L) ) Z? s jf 



K E T (L)-2rE K (L)- ^E T (L)- — P.EMl) + -P;EMLl ) /. - -P EM DP,./, 
-ip-(L)D2j| = fc. 
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Note that the boundary terms from the first set of integration by parts disappear after 
we conflate t with the group parameters. 
Using E X (L) = and the fact that 

(m — 1 or m— 1 or 

L - E Ec-Wfl— K — *- E Ec- 1 )* 25 ^ 7 "™-' 
m=l i=0 m=l i=0 ° T ™ 

(see page 220 of [6]), we can eliminate A. Thus, 

m-l m-l „ r 

A( S ) = 2rP-(L) - L + E + E E (-1)^^, (37) 

m— 1 2— u m— 1 i— u 

where the constant of integration has been absorbed into the Lagrange multiplier and 
hence we obtain two Euler-Lagrange equations in two unknowns. 

To calculate the conservation laws associated to the invariantized Euler-Lagrange 
equations, we must compute the moving frame Ad(p)" 1 and the vector of invariants 
v(I). To compute the former we proceed as in Section [2~3l we calculate the Adjoint 
representation Ad(g) of SE(3) with respect to its generating infinitesimal vector fields 

v a = d x , v 6 = d y , v c = d z , v Q = yd z - zd y , vp = xd z - zd Xl v 7 = xd y - yd x , 

and evaluate it at the frame (|34|) . which yields 

Adipiz))- 1 

where 





( PFS 







{ DXp FS 


D Pfs D 



PFS = ( X S ^ =^=- ) 

is the Frenet-Serret frame, D is the following diagonal matrix 

/ 1 
13= -1 
\ 1 

and X is the matrix 

( ° 

X = z 

\ -y 

As seen in Section [2~3l to obtain the vector of invariants we must first find the 
boundary terms that are linear in the I 2 k- To do so, consider the boundary terms 
obtained from the calculation of the invariantized Euler-Lagrange equations 
(2rE T (L) - A) Jf + (Zf E T (L) - 2ZV S E T (L) ~ V a E K (L)) Jf + (E"(L) + ^E T (L)) V S I% 
+ (kE t (L) - 2rE K (L) - ^E T (L) - ^V s E r (L) + lZ> s 2 E T (L)) J| - \V s E T {L)V s q (38) 
+iE-(i)25?/| = fc. 




Substituting T> S I 2 , T> S I 2 and Cg/f m Q38jl by the differential formulae 

z 2 ) 



v.q = -Ki% + q i + Ti} 



V s P 2 =-Tl y 2 +P 12 , 

V 2 S I* = tkI* - tJI - 2rlf 2 - t<I} + J? 1S 
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which were obtained from ([8]). yields that the boundary terms are linear in the I^ki 

-V S E K (L) - —T> S E T (L) 

K 

E«(L) 



( ij )( - K E«(L)+rE-(L)-A( S ) ) + ( Jf J* ) 



+ ( J% ^12 ^112 ) 



/ is? 



V 2 S E T (L) - —^T> S E T (L) + kE t (Zy) - rE K (L) 



1 



£> S E T (L) 



V 



k. 



Finally, adding the products of the matrices of invariantized infinitesimals Q a (I) with 
the vectors C a yields the vector of invariants 

/ rE-(L)- K E K (L)-A(s) \ 



v(I) 



tE t (L) — kE k (L) — A(s) 
-V S E K {L) - 1V S E T {L) 

E7{L) 



rE K (L) 



V 



E"(L) 



where the f2 Q (/) are 

/ 1 \ 







V o y 

Thus, the conservation laws are 



n x (i) 



( 


\ 


1 























V 


1 / 



n*(j) = 



( 





\ 











1 














K 





1 













PFS 





DXpps 


Dp FS D t 



I tE t (L) — kE k (L) — A(s) 

-V S E K {L) - ^V S E T {L) 
l©2 E r (L) _ ^ x, sE r (L) + kE t (L) 

P-(L) 
-±V B E T (L) 
V E*(L) 



rE«(L) 



Cl 
C2 



c,(39) 



/ 



where Ci = (ci,C2,cs) T , C2 = (C4, c^, cq) t are constant vectors and A(s) is equal to 

(EH). 

We shall see in the remainder of this section how the conservation laws (|39[) can 
help reduce the integration problem. 

To demonstrate this, we will start by simplifying the conservation laws (|39l) in two 
steps. These simplifications will then lead to an overdctcrmined system of equations 
for x, y and z, which will be solved with relative ease once k and r are known. Finally, 
we will give a reason behind the choice in the order in which we solve the equations. 

In the first step of the simplification, we apply an element of SE(Z), say Ad(g) , 
to both sides of Ad(p(z))~ 1 v(I) = c such that it maps ci and C2 to the z-axis. Let 
Ad(g) act on c as follows 

R \f ci A / ci ' 
DAR DRD \ c 2 I cl 



Ad(g)c 
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where R is the three-dimensional rotation 

(cos j3 cos 7 — sin a sin f3 cos 7 — cos a sin 7 — cos a sin /3 cos 7 + sin a sin 7 
cos /3 sin 7 — sin a sin /3 sin 7 + cos a cos 7 — cos a sin /3 sin 7 — sin a cos 7 
sin /3 sin a cos /3 cos a cos /3 

D is the diagonal matrix, D — diag(l, —1, 1), and A the matrix, 

/ -c b 
A= c -a 
\ -6 a 

We can easily verify that the Adjoint representation of SE(3) does not act freely on 
the constant vector c, since it preserves the length of Ci and the quantity cfDc2- 
Indeed to prove the latter, we multiply through 

DARci + DRDc x = c 2 

by cf R T D, and then obtain that 

Ci Aci +c{ DC2 = Ci DC2- 

=0 

Thus, let Ad(g)^ 1 send c to 

c = ( |ci| ^gp ) T . 

Applying Ad(g)^ 1 to the conservation laws yields 

M{g)~ x M{p{*))~ x v{I) = Adig^c, 
which reduces to 

Ad{p{~z))- l v{I) = C (40) 

by the equivariance of the right moving frame Ad(p(z))~ 1 . 

The second step of our simplification consists of applying Ad(p(z)) to (|40"|) to 
obtain the following system of equations 

\ Cl \z s =v^(I), (41) 
l -^z~=v^(I), (42) 



l -^{x s y ss - y s x ss ) = i> (3) (/), (43) 

K 

\cx\(xys ~ yx s ) + ^j^-Zs = v^(I), (44) 

l C l| 

— (x^ s y- yT s x) - Cl , DC , 2 z7s = v^\l), (45) 

K K|Cl| 

— (x(z s x7s - x s z^ s ) - y{y s z7s - zlyTs)) + Cl . DC . 2 (x s y7 s - y s x7 s ) = v^ e \l), (46) 

K K\Cl\ 

where we have used (I) to denote the j-th component of u(J). This overdetermined 
system of equations can now be solved more easily. 
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The two first integrals of the Euler-Lagrange equations are obtained as follows. 
Define B = diag( 1 1 1 ), which satisfies B = Ad(p)~ T BAd(p)~ 1 . Then 
we obtain the first integral of the Euler-Lagrange equations, 

v T (I)Bv(I) = C T BC, 

which is equaivalent to 
(tE t (L) - kE"(L) - A(s)) 2 + (-V S E*(L) - ^V s E T {L)f 

(47) 

+ (iV^(L) - %V S ET(L) + kE t (L) rE-(L)) 2 = c? + c| + c\. 
For the second first integral, we define 

D=( ° J V £ = diag( 1-11), 

which satifies D = Ad(p)~ T DAd(p)~ 1 . The first integral of the Euler-Lagrange 
equations is then 

v T (I)Dv(I) = C T DC, 

i.e. 

(rE-(L) - kE-(L) - X(s)) E T (£) + (-V S E K (L) - ^V S E T {L)) \V S E T {L) 

(48) 

+ {\V 2 S E T {L) - ^2? S E-(L) + kE t (L) - rE K (L)) E K (L) = Cl c 4 - c 2 c 5 + c 3 c 6 . 

We can use the above first integrals to determine n and r. Once we have solved 
for these, we use the system of simplified conservation laws to solve for x, y and z. 
Hence, solving Equation (|41|) gives 



z(s) = j^-r [v^(I)ds. (49) 



|ci| 

Next, multiplying Equation (|43|) by — Cl ^ c ®± 2 and adding it to Equation (|46|) yields 

—(x(z s x7 s -x s z7s) - y(y s z7 s - z s y7s)) = v (6 \i) - Cl , D ^ 2 v ( - 3 \i), 

K |Ci|^ 

which simplifies to 



| Cl | \? B (^(x ■ x) - 1) - -z~ s V s (k ■ x)j = «t;W(J) - ^^(f), 

where ~2?g(x • x) — 1 = x • x^g and ^2? s (x • x) = x • x^. Setting 2? s (x • x) = h(s) and 
substituting £, by (|41[) and by its derivative yields a linear equation for h 



Solving for ft, we obtain 

Ks) = / ^ (2 K (>(/) - ^^ (3) (/)) /^V) + 2) d. 

Hence, 

^(x.x)=^)(I)/-A jy (2 K 
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To solve Equations and (|50)) . we use the cylindrical coordinates 

x(s) — r(s) cos 0(s), y(s) = r(s) sin#(s), z(s) — z(s). 
Starting with Equation (I50p . we obtain 
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r(sf 



as 2> a (x-x) = V s (r{s) 2 + z(s) 
(|44l) becomes 

1 



/i(s)ds — z(s) , 

2\ 



(51) 



r(sre s = — [v^(I)- 



After applying the change of coordinates, Equation 
Ci T L>c 2 



c i 



and then solving for 9 yields 



9(s) 



1 



r( S ) 2 | Cl | 



« w (i)-3^ (1) (J))<i*- 



(52) 



To recover x, y and z, we act on x, y and z as follows 

x4x=Rx|a, 

where R is the three-dimensional rotation 

(cos (3 cos 7 — sin a sin /3 cos 7 — cos a sin 7 — cos a sin /3 cos 7 + sin a sin 7 
cos j3 sin 7 — sin a sin /3 sin 7 + cos a cos 7 — cos a sin /? sin 7 — sin a cos 7 
sin /3 sin a cos /3 cos a cos /3 

and a = ( a b c ) T is the translation vector, with 



tan 



d C 5 Ci 

— c H 



■\/|ci| 2 cos 2 (3 — 
C3 

C2C! T DC2 



7 = tan" 



, c 2 C 4 Ci 
= — c H 



-1 f c 2 c 3 sin/3 -I- Cia/|ci| 2 cos 2 y3 — c§ 
^€103 sin/3 - c 2 v / l c i| 2 cos 2 /3 - c§ 

2 - c 1 c 1 T L»c 2 



C3 



C3 Ci 



C3 C3|ci|" 
and where /3 and c are free. 

Although only four of the equations of the system were used to solve for x, y 
and z, we know that the remaining two equations have been satisfied. Indeed, if we 
differentiate Ad(p(~z))~ 1 v(I) = C with respect to s and multiply by Ad(p(z)), then 
we get 

vMi) = v s (AdipCzmMpwr 1 ^!), 

which is equivalent to 



V a v(I) 








K 














\ 




— K 





T 


















— r 






























— AC 















-1 


K 





— T 




I 





-1 








r 





/ 



V(I). 



(53) 



The above system of equations not only forms part of an elimination ideal as it only 
involves invariants, but because the invariants appear only on the right-hand sides of 
the Equations gl]), g2]), g3]), (|44j) , (g5]) and (01]), they also encode the relationships 
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between the equations themselves. Thus we see V s (Equation (HH)) = ^(Equation (|42p) 
and so forth. Using the Equations in we can eliminate (|4"2")l and from the 
system. 

In the following two well-known examples, we can verify that we obtain the 
expected solutions. 

Example 4.1. For the Lagrangian with L = 1, Equations i51\) and i5ty) yield 

the equation of a line in parametric form as the solution that minimizes the arc length, 
as expected. Note that |x s | = 1 imposes a condition on the constants of integration. 

Example 4.2. Consider the Lagrangian J K 2 ds with torsion, t, equal to zero. Then 
we obtain that k satisfies the Euler- Lagrange equation 

K ss + -K i = 0, 

which is the same equation as for the SE(2) case, or the first integral 

k a + Ak 2 s = |c x |. 
From Equation we know that 



c\Dc 2 = 0, 

and from Equations \51)) and \5ty) we obtain 



z(s) = --, — r / n As, 
c i 



r{sf 



rds 



ds 



k ds 



0{s) =A, 

where A is a constant. So the solution x lies on a plane that includes the z-axis, as 
expected. 



5. Conclusion 



Noether's First Theorem is a well-known result which provides conservation laws 
for Lie group invariant variational problems. In recent work [2], the mathematical 
structure of both the Euler-Lagrange system and the set of conservation laws was 
given in terms of the differential invariants of the group action and a moving frame. 
It is the knowledge of this structure that allows one to solve the invariant variational 
problems under some group action with relative ease. In this paper, we examine 
one-dimensional variational problems that are invariant under the actions of SE(2) 
and SE(3). For both cases, we obtain the invariantized Euler-Lagrange equations 
and their associated conservation laws in the new format, from which we can then 
obtain the solution to the variational problem by quadratures. One can verify that 
this method leads to a far simpler computational problem than the one given in the 
original variables. 
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